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Abstract. Recently, PortcUi ct al (2003) have semi- numerically obtained a 
functional form of the probability distribution of fluctuations in the total energy flow 
in a model for fluid turbulence. This follows earlier work suggesting that fluctuations 
in the total magnetization in the 2D X-Y model for a ferromagnet also follow this 
distribution. Here, starting from the scaling anzatz that is the basis of the turbulence 
model we analytically derive the functional form of this distribution and find its single 
control parameter that depends upon the scaling exponents and system size of the 
model. Our analysis allows us to identify this explicitly with that of the X-Y model, 
and suggest a possible generalization. 



1. Introduction 

Scaling is an important feature of natural phenomena, arising in many degree of freedom 
systems that are highly correlated. In reality these systems support a finite range of 
scales, from the microscopic to the system size. If the number of degrees of freedom is 
sufficiently large, these systems will still fall into the framework of critical phenomena pp. 
Such "inertial" [2j systems include a disparate range of phenomena and, non-intuitively, 
have recently been suggested to have a common signature in the statistics of fluctuations 
in global measures of activity. This probability distribution function (PDF) has been 
compared numerically [2] for a range of models including some for out of equilibrium 
critical phenomena, notably a sandpile, a forest fire model, a depinning model and a 
stacking model for granular media. After normalization to the first two moments these 
PDF were found to collapse onto that of fiuctuations in models for equilibrium critical 
phenomena. 

The functional form of this curve has been identified semi-numerically jS] with the 
distribution [2j: 



where the constants iT, h and s are fixed by setting the moments Mq = 1, Mi = and 
M2 = 1, leaving a single parameter, ag. 

The statistics of fiuctuations in global quantities have been explored both 
experimentally and theoretically for fiuid turbulence in closed systems. In the 
experiments of Labbe et al. ^ Ej, the normalized PDF of fiuctuations in the power 
provided to both rotor blades stirring a closed cylinder of gas at constant angular 
frequency was also found to have collapse onto over a range of Reynolds number 
Rg. These results have been compared with the PDF of fiuctuations in the total 
magnetization in the 2D X-Y model for a ferromagnet [Hj which also has been identified 
with f3]. In the experiment reported by [Zj (their Figure 2), the normalized 
PDF of fiuctuations in wall pressure, rather than injected power, appear to follow 
these non- Gaussian statistics with insensitivity to Reynolds number [7.. Recently, 
jH] treated a model for closed turbulence and obtained a family of curves of the form 



PACS numbers: 05.40.-a,05.50-|-q,47.27.Eq,64.60.Ak 
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(0) semi- numerically for global fluctuations of kinetic energy which within a sign the 
experimentally measured total power is believed to act as a good proxy [8]. 

There is a continuing debate (see [0111111111111211131111113112] and references therein) 
concerning the origin of this apparent universality 12] • The aim of this paper is not 
to establish the existence, or lack thereof, of a universality class. Furthermore, we do 
not address the correspondence between experimentally measured quantities and those 
captured by models for turbulence (see 0). 

Here, we analytically derive for the model for intermittent turbulence in a flnite 
sized system treated by [8j. We obtain function of the model parameters. The 

analysis then leads to a direct identiflcation with results obtained previously for the 
2D X-Y model ^Tj, elucidating the origin of the value ~ tt/2 obtained for that 
system jH] E] • We then suggest that the features of the model that are intrinsic to this 
calculation are rather generic and discuss how they may encompass the wide variety 
of systems which have also been previously identifled as exhibiting the same functional 
form for the fluctuation PDF [2]. 

2. Model for turbulence in a finite sized system 

Portelli et al. [8^ obtained semi-numerically for intermittent turbulence in the 
framework of the K062 hypothesis which models fluctuations in the energy in the flow. 

The essential features of this model are structures on a range of length scales 
from a smallest size /i = r/ to the system size 1^ = L, corresponding to the 
dissipation and driving length scales respectively. The Reynolds number of the flow 
is then deflned as Re = (L/r])^^^ [T^, and the ratio between successive lengthscales 
(Ij/lj-i) = so that = {L/r]f. 

Following jH] we wish to calculate the statistics of the total energy in the flow: 



from a model expressed in terms of an intermittent energy transfer rate ej on lengthscale 
Ij drawn from a PDF which has moments: 



e{t) = ei(t) + e2(t) + ■ ■ • + e,-(t) ■ ■ ■ + e^(t) 



(2) 




(3) 



with the condition < ej >= eo which flxes /i(l) = 0. It follows that the standard 
deviation: 




The individual are assumed independent, each with PDF Pj{ej), giving: 




(5) 
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where Pi{k) = J Pi{ei)deiexp{-ikei). 

In the framework of K062 this scahng system supports a cascade from large to small 
scales, with the intermittency parameter r(2) = — /i(2). Importantly, although we can 
envisage a cascade, the above model does not explicitly require one and will map onto 
other models provided that the basic assumptions, namely of ej that are independent 
and drawn from a PDF with the scaling property ()3I4|) . hold. 

The Pj will depend upon the details of the system, to make progress we first consider 
a tractable choice, the distribution: 

P,(6,) = A,e;^~'e-^^^^ (6) 

and later explore how this may be generalized. The influence of the microscopic 
distribution Pj has been explored in the context of the 2D X- Y model in jTHj. This 
choice for P{ej) was used to evaluate (0) semi numerically in j8|; we will now evaluate 
it analytically and as a corollary directly obtain the solution previously obtained semi- 
numerically for the 2D XY model jS] . 



3. Evaluating the Characteristic Function 

The constants Aj, uj and aj of © are fixed through Q so that Eq. (jSJ can be written 
as 0: 

1 

Pie) = — / e'^'dke-^'' where (7) 



277 J-c 

7 = 1 J J 7 = 1 



Af 1 TV 1 n 



where 



/i = = exp(ja) - 1 (9) 
d = (/i(2)/3) In A and where we define total variance 



N N 



^- = E^J = ^oE/. (10) 



- E = ^0 

i=i i=i 

and P = Neo/a. For a specific system, if one has the values of yu(2) and A, and that (0) 
is a good approximation for Pj{ej), one can evaluate (jHI) numerically |H]. Here, however, 
we wish to establish why appears to also describe the 2D X-Y model. We proceed 
by analytically evaluating (jH}, and begin with Sn- By formal expansion: 

We need to find F„ = J2f=i fj" with the condition that Fq = N and Fi = J2f=i fj = 
{N/jS)"^. From our definition of the fj we can evaluate Fi and obtain 

exp(A^a) = 1 + (A^ + Ny(3^){l - exp(-a)) (12) 
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In the limit of N / (3"^ > 1, that is, ^ 1, so that | a |<^ 1 for finite system size L (ie 
A ^ 1) this gives: 

e^" ^ 1 + ^ (13) 
We can now expand the F„ in e^"- and substitute (fT!?jl : 

N N N 

F„ = ^(e^" - 1)" = ^ e"^" - n ^ g^""^)^'^ + . . . (14) 
j=i j=i j=i 

^nNa g(n— l)7Va 

— n 7 -— + . . . 



Q~na ]^ g— (n— l)a 

and in our hmit ^ 1, | a |-C 1 such that > 1 we can take | a | //?^ ~ 1 to 

give to lowest order ~ {N'^ / {na)){{Na) / jS'^Y where n>l {Fq = N). Using this in 
(fTT|l gives 

^ = iap -iaZ\n(l + ika/Z) (15) 
dk 

where Z = j3/{Na). This now readily integrates to give: 

SNik) = ika[3 + ipik) where (16) 
9 / cr \ , / .a 



i\){k) = ikaZ -Z'i^l + tk—j In (^1 + ik—j (17) 

With (j) = e — e^N = e — Pa we can write (jH} as: 

1 

P{e) = — / e""f'e-^^''^dk (18) 

The limit Na — > 0, 2 — > oo, corresponds to retaining terms up to k"^ in ()17|) and 
immediately gives a Gaussian distribution for ^(e). It is tempting to take this as the 
Gaussian limit of P{b), a — 0, or yu(2) 0. However, since we insisted that the PDF 
Pji^j) are scaling, this limit would yield (Xj — > in The case where the Pji^j) all 
have the same (that is, non-scaling) aj, corresponding to Na ^ as this gives F2,,n 
above. To evaluate (fTH|) in general we need to retain the property of scaling cxj, thus 
excluding this limit. 

We evaluate (fT^ by the method of steepest descent: 

P{e) = ^e-^(^+^)-^'^"^ (19) 
crv27r 



which is of the form with ag = Z'^ + 1/2, that is, 

and u = — e(a/eo) + {Nd + A) with A just given by the normalization constant K. 
The present choice of Pj{ej), Eq. ©, gives exactly: 

^^^''^ = (21) 
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with the Kj and 7^- fixed by the moments. Other forms of Pj{k) for which (|2H) is a good 
(Pade type) approximant will yield a P{e) of the form (^. We write ()21|) as: 

00 (-ik)P 

^#) = E^<^?> (22) 
p=o 

~l+zA;eo-— ^(1 + -) + ••• 
2 7i 

given that < ej >= eo = i^jlj- All the coefficients in this expansion are fixed if we insist 
that any Pj that we consider has the same scaling (jlj) for aj, so that exp(ja) = (1 + 1/7^). 
For the PDF (jH)) this gives 7^ = uj which is exact. For example, one might anticipate 
that in a correlated system that local fiuctuations may be multiplicative. An appropriate 
model for multiplicative noise is a lognormal PDF: 



Pj{ej) = —=- — exp 



2-1 



(23) 



which has Pade type approximant of form (|2H) fixed by: 

< >= tPe\p'-^\ 1 + - = e"^^ (24) 

Thus if this is a good approximant, the lognormal also yields a curve of form (|H). 
4. Results for the turbulence model. 

To make a direct comparison with the results of [HI IHI we write ()20|) in terms of the 
parameters relevant to the turbulence model: 
1 3 / M2i 

SO that ag depends weakly on the Reynolds number of the fiow, on (experimentally 
determined) r(2) and through the logarithm, on the free parameter A. In Figure 1 we 
show normalized curves for the range Re = [10^, 10^, 10^] corresponding to ag = [7. 7.. 3. 6] 
explored by 0, A = 2 and typical values of r(2) |TH]. Curves for r(2) = —0.25 are 
shown in the main plot and for r(2) = —0.2 (used jHI) in the inset. In both cases these 
show good correspondence with the solutions from the model (see Figures 1 and 2 of [S]). 
The curves also fall close to each other explaining the relative insensitivity to Reynolds 
number jH], and also fall close to that for ag = n/2 identified by for the X-Y model, 
which is also plotted and which we discuss next. 

A more sensitive method for identifying ag from data is to calculate the third 
moment of (0) directly [201 rather than compare curves. With Mq = 1, Mi = 0, M2 = 1: 

Na I rlnA \ / 



SO its magnitude increases with decreasing ag which from ()25|) corresponds to increasing 
Re, consistent with data. The handedness is not determined here; this depends upon 
how the global quantity measured experimentally relates to that of the model; namely 
the energy fiow within the fiuid 0. 
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Figure 1. Function |(TJ is plotted for Og = [4.24,2.76,1.90] (blue lines), Og = 1.74 
(red line) and Ug = tt/2 (diamonds) relevant to the model 8 for a closed turbulence 
experiment, the 2D X-Y model and [2] respectively. The inset shows the same axis 
ranges and replaces the blue curves with Og = [7.76,5.18,3.62] relevant to IH|. The 
handedness is chosen for comparison with [21 El |H1 • 



5. The 2D X-Y Modeh 

The total magnetization in the spin-wave approximation to the 2D XY model was also 
shown for a range of system sizes to collapse onto ((T)), with ~ 7r/2 (for a study of 
the full 2D XY model see |2I])- The PDF of total magnetization M can be written as 

Q{M) = r JLe-^^'^e^ (27) 



' — oo 

oo , 



2'Ka 



2V a„ 1^1 



where 5 = ^ -^^\/ " TT' S = TTp " 

and 7q specifies the lattice Green's function. Bramwell et al jH] evaluated this numerically 
to demonstrate that it is well described by ((H). We now show that the sum S is related 
to the sum Sjy (jHl) by writing: 

^ ln(l + zkeoU) ^ 1 ^ Hkeof.r ...^ 
-^N = -}_^ 7 = -r 2^ (28) 

p=l JP p=l JP m=l 

where < e >= Ne^. If we then make the identification: 

p 



ar. 



2 V 



then Q{M) (P7|) has the same functional form as P{s) so that it also shares the 
same distribution ^ to within the approximations made here, namely, that following 
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expansion (14) we have neglected terms of order and that (27) is also an 

approximation, good for large. Importantly, from the definition of in ()27|) . the 
r.h.s. of ()29|1 is independent of Ng. In this sense we have approximately evaluated the 
integral (jTfj) in the thermodynamic limit. 

It then just remains to estimate ag for the 2D X-Y model. In (3, this was 
approximated asymptotically, here we simply note that insisting that the normalized 
Q{M) and P{s) share the first three moments yields ag = 1/2 + {a2/2)^/al from (j^Bjl 
and equation (21) of [3 . They also calculated the normalized third moment for a square 
lattice. Their value of M3 = —0.8907 gives ag = 1.7428 which will give curves close to 
those for ag = 7r/2 as shown on Figure 1. Our analysis is thus consistent with both a 
value of ap ^ 7r/2 0, and the asymptotic exponent of |3|. 

6. Generalization 

A variety of disparate systems have recently been shown numerically 12] to have a 
common signature in the statistics of fluctuations in a global measure of activity which is 
of the form (^J). These include out of equilibrium critical phenomena, notably a sandpile, 
a forest fire model, a depinning model and a stacking model for granular media. 

We will now argue that the scaling anzatz which was our starting point for the 
model for fluid turbulence in section 2 and our derivation of may also encompass 
these disparate systems. 

The ansatz we chose corresponds to that of a scaling system that generates spatial 
structures or domains (patches) on length scales from a smallest size h = t] to 

the system size 1^ = L. "Length scale" in this more general sense means "appropriate 
characteristic measure" i.e. length in one dimension, area in two dimensions or volume 
in three dimensions. In a dynamical out of equilibrium system, such as a sandpile or a 
forest fire model, a steady state is achieved by driving on the smallest length scale h = rj 
and by means of open boundaries, removing structures on the system size L. In section 
2 we considered a system driven on the largest scale L and dissipating on the smallest, 
mapping onto fluid turbulence in a closed system. In a model for a ferromagnet, the 
system may fluctuate about an equilibrium, but nevertheless has a minimum patch size 
(one spin), a maximum patch size (the system size), and scaling of patches in between. 

The global quantity e is now taken to be associated with the total number of 
instantaneously active sites within each patch. In a model realized numerically, such as 
a forest fire or avalanche model, instantaneously active sites are those seen at a given 
timestep in the computation. In a forest fire model, active sites correspond to burning 
trees, in an avalanche model, to relaxing sites in evolving avalanches j2]. The global 
quantity may refer to the energy dissipated by these sites, or simply refer to the time 
evolution of their spatial distribution as in the case of space occupied by anisotropic 
particles settling under gravity or magnetization of spins in a ferromagnet [2j. 

The common feature of these systems is that at any instant in time there will be 
rrij (t) patches on any length scale Ij and associated with each patch, e* of this quantity. 
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On each length scale Ij we then have = mj(t)e* and in total, e given by 0- We 
now assert that in common with the turbulence model, the e* are independent and have 
intermittent, scaling statistics (jH)). We can then envisage the following generic scaling 
system which comprises: 

(I) Non space filling, intermittent patches: The details of < > depend on the 
system, for example the probability of patches Zj_i merging, and/or patches Ij^i breaking 
up to form patches on Ij. We take as a necessary condition of scaling that the moments 
obey: 

< m] > if'^ =< m]„i > if I =< m% > (30) 

If the system were space filling, 7(1) would be 1 so 7(1) < 1 implies non-space filling 
patches. Allowing 7 = 7(5') permits intermittency. 

(II) Fractal support: On any patch there will be a density of active sites e*/lj which 
in general can vary with Ij; for a system which is scaling we can however take: 

A = = f31) 
If lf_i ^ ^ 

where a = 1 is the special case of uniform density on all patches, and patches that do 
not have fractal boundaries. 

(III) Conservation: Scaling implies that there is no preferred Ij on which the active 
sites accumulate so that the mean will be just the ensemble average determined on any 
length scale. This is consistent with conservation of active sites when patches merge 
{Ij Ij+i) or break up {Ij — > Ij^i). 



It follows from (I), (II) and (III): 




(ag-7(g)) /, ^ -/i{q) 

(32) 




The condition < ej >= to fixes 7(1) = a or yu(l) = 0. The details of the system 
specify fj,{2) which then fixes the standard deviation of © expressed through (jH^ and 
immediately leads to (jH). 

(IV) Finite size: We finally specify the number of length scales A^; given scaling, a 
choice is constant {Ij/lj^i) = so that = {L/r])^. 

In summary then, this scaling ansatz is that < tj >= cq, < e| >= Cq (/j/L)^^*-^-* 
with /i(2) 7^ and {Ij/lj-iY = ^- ^"^Y system that is specified by this anzatz and is 
well approximated by (pT|l will share the same behavior in the statistics of global 
activity P{e) that we have calculated above for the turbulence model. Importantly, 
these conditions may apply to more than one quantity in a given system, and any such 
quantity will share these same statistics. 

For a given system, the curve (P) is specified by which is a function of the 
system parameters N ^ /i(2) and A. This family of curves is however insensitive to Og 
|2Uj . This, combined with the practical difficulty of obtaining good statistical resolution 
over fluctuations ranging over several orders of magnitude suggests a straightforward 
reason for the close, but not exact, curve collapse that has been reported in figure 2 of 
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[2]. Importantly, we do not extend this argument to the 2D X-Y model; rather in this 
case we have utilized the correspondence of (7) with the result of [17] (equation (27)). 

7. Summary. 

From the starting point of a model for fluid turbulence in a finite sized system, previously 
treated semi- numerically by [S] , we have analytically derived the functional form of the 
PDF of global energy flow in the system. This yielded the dependence of its single 
control parameter Qg on the intermittency parameter, the ratio between lengthscales, 
and the smallest and largest scale lengths in the system (i.e. the Reynolds number). 
We then directly identified this function with that previously obtained for fluctuations 
in total magnetization in the 2D X-Y model and thus elucidated the origin of the 
previously identified value ~ 7r/2 |i2j. The PDF was shown to be relatively insensitive 
to variations in a^, explaining the previously reported close correspondence of these 
curves for the turbulence model and the 2D X-Y model E] • 

Importantly, the functional form of the PDF that we derive is just that of the 
sums of a large but finite sets of independent numbers drawn from PDF with moments 
that are scaling. This corresponds to a model of intermittent turbulence in which one 
also envisages a cascade, but the cascade property is not intrinsic to the calculation. 
We suggest that this system is rather generic and may encompass the wide variety of 
systems which have also been previously identified as exhibiting the same functional 
form for the fluctuation PDF 0. 
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